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Abstract. We describe results from semi-analytical modelling of star formation in protocluster clumps of 
different metallicities. In this model, gravitationally bound cores form uniformly in the clump following a 
prescribed core formation efficiency per unit time. After a contraction timescale which is equal to a few times 
their free-fall times, the cores collapse into stars and populate the IMF. Feedback from the newly formed OB 
stars is taken into account in the form of stellar winds. When the ratio of the effective energy of the winds 
to the gravitational energy of the system reaches unity, gas is removed from the clump and core and star 
formation are quenched. The power of the radiation driven winds has a strong dependence on metallicity and 
it increases with increasing metallicity. Thus, winds from stars in the high metallicity models lead to a rapid 
evacuation of the gas from the protocluster clump and to a reduced star formation efficiency, as compared 
to their low metallicity counterparts. We derive the metallicity dependent star formation efficiency per unit 
time in this model as a function of the gas surface density E g . This is combined with the molecular gas 
fraction in order to derive the dependence of the surface density of star formation Y^sfr on E s . This feedback 
regulated model of star formation reproduces very well the observed star formation laws in galaxies extending 
from low gas surface densities up to the starburst regime. Furthermore, the results show a dependence of 
EsFij on metallicity over the entire range of gas surface densities, and can also explain part of the scatter 
in the observations. 
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1 Introduction 

Over the last two decades, the dependence of the star formation rate surface density (Tisfr) on the gas surface 
density (S 9 ) and eventually on other physical quantities has been extensively investigated both observationally 
(e.g., Kennicutt 1998; Wong & Blitz 2002; Boissier et al. 2003; Bigiel et al. 2008; Blanc et al. 2009; Onodera et 
al. 2010; Tabatabaei & Berkhuijsen 2010; Heiner et al. 2010; Schruba et al. 2011; Bolatto et al. 2011) as well 
as theoretically and numerically (e.g., Tutukov 2006; Krumholz & Thompson 2007; Fuchs et al. 2009; Silk & 
Norman 2009; Krumholz et al. 2009a; Papadopoulos & Pelupessy 2010; Gnedin & Kravtsov 2011; Narayanan 
et al. 2011; Fcldmann et al. 2011; Vollmer & Leroy 2011; Braun & Schmidt 2011; Monaco et al. 2011; Kim 
el al. 2011; Dib 2011a, b). Determining the rate of star formation in a given tracer of the gas surface density 
requires quantifying the fraction of that tracer as a function of the global gas surface density and a description 
of the efficiency at which the star forming gas is converted into stars per unit time. The relationship between 
^SFR and the surface density of the molecular hydrogen gas Y*h 2 is given by: 

SsFi?. = S g Jh 2 -, (1.1) 

T 

where fn 2 = ^H 2 /'^g is the molecular hydrogen mass fraction is the atomic- molecular star forming complexes, 
and SFE T is the star formation efficiency over the timescale t. Krumholz & McKee (2005) proposed a theory 
in which supersonic turbulence is the dominant agent that regulates star formation in giant molecular clouds 
(GMCs). They derived a core formation efficiency per unit free-fall time, CFEff, which is given by CFEff = 
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0.15a-° r 68 A4-°- 32 , where a mr and M are the virial parameter and the rms sonic Mach number of the GMC, 
respectiveljQ. By assuming that only a fraction of the mass of the cores ends up in the stars, this CFEff 
can be converted into a star formation efficiency SFEff = rj x CFEff (rj < 1). An alternative theory has 
been recently proposed by Dib et al. (2011) and Dib (2011a, b) in which the star formation rate in protocluster 
clumps is primarily regulated by feedback from massive stars and in particular through energy injection in the 
clumps by stellar winds. In the following sections, we briefly describe the main constituents of this model. 

2 Feedback regulated star formation 

The model follows the formation of dense gravitationally bound cores in a protocluster clump. Cores form in 
the clump with a given core formation efficiency per unit time and follow a local mass distribution that is the 
result of the gravo-turbulent fragmentation of the clump. In their series of models, Dib et al. (2011) varied the 
core formation efficiency per unit free-fall time (CFEff) between 0.1 and 0.3. This is consistent with the range 
of CFEff measured in numerical simulations which describe the gravo-turbulent fragmentation of magnetised, 
turbulent, and self-gravitating molecular clouds (e.g., Dib et al. 2008; Dib et al. 2010a). The gravitationally 
bound cores that are generated at every epoch have a mass distribution that is given by the gravo-turbulent 
fragmentation model of Padoan & Nordlund (2002). In this work, we leave out, for simplicity, the role played 
by gas accretion and coalescence in modifying the mass distribution of the cores. The interested reader is 
referred to Dib et al. (2007) and Dib et al. (2010b) for such models. Cores contract over a lifetime which is 
a few times their free-fall time before collapsing to form stars. Feedback from the most massive stars (M* > 5 
M Q ) is taken into account in the form of stellar winds. The formation of cores in the protocluster clump, and 
consequently star formation, are terminated whenever the fraction of the wind energy stored into motions that 
oppose gravity exceeds the gravitational energy of the clump. In order to calculate reliable estimates of the 
feedback generated by metallicity dependent stellar winds, we proceed in two steps. In the first step, we use 
a modified version of the stellar evolution code CESAM (see appendix 1 in Piau et al. 2011) to calculate a 
grid of main sequence stellar models for stars in the mass range [5-80] M Q (with steps of 5 M©) at various 
metallicities Z/Zq = [1/10, 1/6, 1/3, 1/2, 1, 2] (Z Q — 0.0138). The evolution of massive stars is followed using 
the CESAM code for ~ 1 Myr, on the main sequence. The characteristic stellar properties, which are the 
effective temperature T e //, the luminosity L*, and the stellar radius i?* are then used in the stellar atmosphere 
model of Vink et al. (2001) in order to calculate the stellar mass loss rate M+. Vink et al. (2001) did not derive 
the values of the terminal velocities of the winds (foo), therefore, we use instead the derivations of Woo obtained 
by Leitherer et al. (1992). 

The power of the stellar winds is given by M+v^ . This quantity is displayed in Fig. [1] for the models with 
different metallicities. The values of M+v^ are fitted with fourth order polynomials (overplotted to the data) 
and whose coefficients are provided in Dib et al. (2011). The M+v^ — M+ relations displayed in Fig. [Hallow for 
the calculation of the total wind energy deposited by stellar winds. The total energy from the winds is given 
by: 

s~ . r r i2 ° M " ( w(M - )A y M - ),, M *• ("> 

Jt'-o iif,=5 m y * J 

We assume that only a fraction of E W i nc i will be transformed into systemic motions that will oppose gravity 
and participate in the evacuation of the bulk of the gas from the proto-clustcr clump. The effective kinetic wind 
energy is thus given by: 

Ek,wind — K E w i nc [, (2.2) 

where n is a quantity < 1 (in this work, we use k = 0.1 for all models). Ek^wind is compared at every timestep 
to the absolute value of the gravitational energy, E grav , which is calculated as being: 

16 f R " 

E g rav = -— tt 2 G / p c {rfr A dr. (2.3) 



3 



o 



*Padoan & Nordlund (2011) found a different dependence of the SFR on a v i r and Ai. The results of their numerical simulations 
suggest that the SFR decreases with increasing otm T but also that it increases with increasing M. 
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Fig. 1. The power of the stellar winds, or wind luminosities, for stars in the mass range 5-80 Mq on the main sequence, 
and for various metallicities. The stellar mass loss rates have been calculated using the stellar characteristics (effective 
temperature, stellar luminosity and radius) computed using the stellar evolution code CESAM coupled to the stellar 
atmosphere model of Vink et al. (200f ). The values of Voo have been calculated using the derivation by Leitherer et al. 
(f992). Over-plotted to the data are fourth order polynomials. The parameters of the fit functions can be found in Dib 
et al. (2011). Adapted from Dib et al. (2011). 



Since higher metallicity stellar winds deposit larger amounts of energy in the clump than their lower metallic- 
ity counterparts, this leads them to evacuate the gas from the clump on shorter timescales. This in turn quenches 
the process of core and star formation earlier and sets a smaller final star formation efficiency, SFE exp . Fig. [2] 
displays the dependence of SFE exp and of the expulsion time, t exp (expressed in units of the free-fall time £//), 
as a function of metallicity for clumps of various masses. 



3 Metallicity dependent star formation laws 

Using the above described model, it is possible to derive the dependence of Tisfr on E s . The star formation 
rate surface density in the feedback regulated mode of star formation is given by: 



(SFE eX p) 
{^exp) 



^SFB, — LgJH2 ~ \ , (o.lj 



where (SFE exp ) and (t exp ) are, respectively, the characteristic SFE exp and the epoch at which gas is expelled 
from the protocluster region for the clump mass distribution associated with a given E s . Writing (t exp ) in terms 
of the clump free-fall time (£//), Eq. l3.1l becomes: 
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Fig. 2. Dependence of the quantities SFE exp (final star formation efficiency) and n exp = t exp / t ,, (ratio of the expulsion 
time to the free-fall time) for selected values of the protocluster forming clump masses and metallicities. These results 
are based on the models of Dib et al. (2011). 



(SFE exp ) 1 

"" 2 '^ XP ) (tff) 
(/*■//) 
{tff) 



^Ih^-^T- (3.2) 



where /*,// is the dimensionless star formation efficiency and which corresponds to the mass fraction of the 
molecular gas that is converted into stars per free-fall time tff of the clumps. (/*,//) and (tff) represent 
characteristic values of f*j / and tff for the spectrum of clump masses found in the GMC for a given value of S fl . 
The quantity fn 2 is the mass fraction of the total gas that is in molecular form. In this work, we use the functional 
form of fu 2 obtained by Krumholz et al. (2009b) who derived fn 2 as a function of the gas surface density and 
metallicity (see their paper or Dib 2011 for the detailed formula), (tff) can be approximated by the free-fall 

time of the clump with the characteristic mass tff(M c h ar ) — 8S c( 3 ^ 4 M^ r 6 Myr where M c i iar fi — M c har/10 6 
M©. The characteristic mass M c h ar is given by : 
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Fig. 3. Star formation efficiency per unit free-fall time in the protocluster clump in the metallicity-dependent feedback 
model. The star formation efficiencies per free-fell use a core-to-star efficiency conversion factor of 1/3. The left panel 
displays /*,// as a function of both M c i and Z — Z/Zq in the original data, while the right panel displays the analytical 
fit function to this data set given in Eg. 13.51 Adapted from Dib (2011a). 



Mr 



max(M < .i irrla ,z.,M G Mc) 



M cl N(M cl )dM cl , 



(3.3) 



where N{M c i) is the mass function of protocluster forming clumps which we take to be N(M c i) — A c iM cl 2 , 
and A c i is a normalisation coefficient given by A c i f^^ Mci ' maa: ' MGMC '' N(M c i)dM c i — e, where < e < 1 is 
the mass fraction of the GMCs that is in protocluster clumps at any given time. In this work we use e = 0.5. 
The minimum clump mass Mri, m in is taken to be 2.5 x 10 3 M© (this guarantees, for final SFEs in the range of 
0.05-0.3 a minimum mass for the stellar cluster of ~ 50 M©) and the maximum clump mass is 10 s M©. The 
characteristic GMC mass is determined by the local Jeans mass and is given by: 



M GMC = 37x 10 6 
Fig. 2] (top) displays M c har a s a function of S, 



(3.4) 



^85 M© pc- 2 , 

, g . The quantity /*,// = SFE exp /n exp is displayed in Fig. [3] 
(left panel) as a function of mass and metallicity (Z — Z/Z®). These models use a value of CFEff = 0.2 
and standard clump and core parameters (see Dib et al. 2001 and Dib 201 la, b for more detail). A fit to the 
(M c i , Z ) data points with a 2-variables second order polynomial yields the following relation shown in Fig. [31 
right panel: 



Uj f {M ch Z') = 11.31 - 4.311og(M cl ) + 0.41[log(M cl )] 2 
-8.28Z + 3.20Z'log(M cl ) - 0.32Z'[log(M cl )] 2 
+2.30Z 2 - 0.89Z' 2 log(M cl ) + 0.08Z' 2 [log(M cl )] 2 . 
Using Eg. 13.51 it is then possible to calculate (/+.//): 

rmax(M c i imax , Ma m c ) 

(/*,//) (Z ,T, g ) = / U Jf (M cl ,Z)N(M c i)dM cl . 



(3.5) 



(3.6) 



Mr. 
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Fig. 4. Characteristic clump mass as a function of the gas surface density (Eq. 13.31 top panel) and the star formation 
efficiency per unit free-fall time in this feedback regulated model of star formation (lower panel). Adapted from Dib 
(2011a). 



Fig. [4] (bottom) displays (/*,//) [Z , £ 9 ) for values of Z in the range [0.1 — 2}. We assume that there is a 
critical value of S g = 85 M Q pc~ 2 below which clumps are pressurised by their internal stellar feedback, such 
that Y>d — Y,g tC rit where E 9 < £ SjC rit and £ c z = Egmc = E 9 when E g > £ ff , C rit- With the above elements, the 
star formation law can be re-written as: 



-sfr - J^fH 2 (Sg,C,Z )T,g 

U*,f. f )(z') 



j 1 / 4 ' 85 M pc 

char, 6 

) ( >:» V i/M . s 



^2 < 1 



./.f)(z') f 



jfi/* \^85 M oP c- 2 J ' 85 M pc- 2 



> 1 



(3.7) 



where T,sfr is in M Q yr 1 kpc 2 , M c h ar is given by Eq. 13.31 and (/*,//) by Eqs. 13.51 and !3. 61 Fig. [5] displays the 
results obtained using Eq. 13.71 for S g values starting from low gas surface densities up to the starburst regime. 
The results are calculated for the metallicity values of Z = [0.1, 0.3, 0.5, 1, 2]. The results are compared to the 
sub-kpc data of Bigiel et al. (2008,2010) and to the normal and starburst galaxies results of Kennicutt (1998). 
The models fits remarkably well the observational results over the entire range of surface densities. Furthermore, 
the segregation by metallicity extends beyond the low surface density regime up to the starburst regime where 
a segregation in metallicity of ~ 0.5 dex is observed. 



The dependence of star formation rates in galaxies on metallicity 



243 




1CT 1 10° 10 1 10 2 10 3 10 4 10 5 

S g (M pc- 2 ) 

Fig. 5. Star formation laws in the feedback-regulated star formation model. Overplotted to the models are the normal 
and starburst galaxies data of Kennicutt (1998) and the combined sub-kpc data (4478 subregions) for 11 nearby galaxies 
from Bigiel et al. (2008,2010). The Bigiel et al. data is shown in the form of a 2D histogram with the colour coding 
corresponding, from the lighter to the darker colours to the 1,5,10,20, and 30 contour levels. The displayed theoretical 
models cover the metallicity range Z — Z/Zq — [0.1, 2]. 

4 Conclusions 

We have presented a model for star formation in protocluster clumps of difference metallicities. The model 
describes the co-evolution of the dense core mass function and of the IMF in the clumps. Cores form uniformly 
over time in the clumps following a prescribed core formation efficiency per unit time. Cores contract over 
timescales which are a few times their free fall time before they collapse to form stars. Feedback from the 
newly formed OB stars (> 5 M Q ) is taken into account and when the ratio of the cumulated effective kinetic 
energy of the winds to the gravitational energy of the system (left over gas+stars) reaches unity, gas is expelled 
from the clump and further core and star formation are quenched. The radiation driven winds of OB stars 
are metallicity dependent. Metal rich OB stars inject larger amount of energy into the clump than their low 
metallictiy counterparts and thus help expel the gas on shorter timescales. This results in reduced final star 
formation efficiencies in metal rich clumps in comparison to their low metallicity counterparts. Both the final 
star formation efficiency and the gas expulsion timescales are combined for a grid of clump models with different 
masses and metallicities in order to calculate the star formation efficiency per unit time (/*,//) in this feedback 
regulated model of star formation. We calculate the characteristic value of /*,// for a clump mass distribution 
associated with a gas surface density, E g . This is combined with a description of the molecular mass fraction as 
a function of S g and the assumption that there is a critical surface gas density (£ g = 85 M pc -2 ) above which 
the protocluster clumps and their parent giant molecular clouds switch from being pressurised from within by 
stellar feedback to being confined by the external interstellar medium pressure. The combination of these three 
elements allows us to construct the star formation laws in galaxies going from low gas surface densities up to 
the starburst regime. Our models exhibit a dependence on metallicity over the entire range of considered gas 
surface densities and fits remarkably well the observational data of Bigiel et al. (2008,2010) and Kennicutt 
(1998). This dependence on metallicity of the KS relation may well explain the scatter (or part of it) that is 
seen in the observationally derived relations. 
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